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1. Introduction 

In the quest for the understanding of the (now not so) new evidence about the increase 
in the expansion rate of the Universe dubbed cosmic acceleration, modifications of the 
structure of Einsteinian gravitation have played a major (and highly debated) role [HE]- 
Although it is still not clear if this unexpected behavior of the Universe should or could 
be ascribed to new aspects of the gravitational interaction, so far a great number of 
modified theories of gravitation have been proposed. The most studied classes of such 
extensions are the ones that are often called "extended theories of gravity". These 
theories are in some sense "minimal" modifications of GR and are designed to be fully 
metric-affine i.e. such that they are completely determined by a metric and a connection 

In this work we will focus on one type of such extensions of General Relativity (GR), 
which takes the name of /(i?)-gravity [1] and it is characterized by an action which is 
non-linear in the Ricci scalar. These theories emerged for the first time in quantum 
field theory in curved spacetime as an effect of the renormalization of the matter stress 
energy tensor [5] and were successively recovered as low energy limit of fundamental 
theories like string theory [6]. Many aspects of the cosmology of /(i?)-gravity have been 
investigated with a variety of techniques, but most of the studies were focused on the 
homogeneous and isotropic case. Even in this simple case the differences between the 
cosmology of Relativity and the ones of these models can be quite dramatic. It is known 
for example that in the some simple models of this type perturbations can grow even in 
a state of accelerated expansion [7j. 

All of these extended models of relativistic gravitation are usually designed 
considering matter as a purely classical system, neglecting quantum properties like, 
for example, the spin. Nonetheless the spin is an important property of the fermions, 
which are the most common type of fundamental constituents of standard matter. It is 
natural then to attempt the introduction of this feature in the mathematical structure 
of relativistic gravitation. In the context of GR it has been proven that one way to 
achieve this goal is to include an additional geometric property of spacetime: torsion 
[H [U [To] . Geometrically torsion is defined as a property of the spacetime that prevents 
infinitesimal parallelograms to be closed. This idea translates directly in the existence of 
an antisymmetric part of the connection and it is the cornerstone of the Einstein-Cartan- 
Sciama-Kibble (ECSK) theory. Since the choice to have only a completely symmetric 
connection was an assumption made in the first years of GR, we can think about ECSK 
theory as a generalization rather than a modification of Einstein gravitation. The 
presence of torsion beside curvature, allows both the spin and the energy of matter 
fields to be coupled with gravitation and opens naturally the way to the construction of 
an Einstein-Dirac theory. 

The inclusion of torsion is therefore justified by the presence of the spin, which is 
itself justified since we know about the existence of fermions. Within fermionic field 
equations in general, and for the Dirac field in particular, the coupling spin-torsion 
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gives rise to fermionic self- interactions [12] . In the framework of ECSK theory, however, 
since the gravitational Lagrangian is the usual Hilbert-Einstein one (with torsion) 
containing only the gravitational Newton constant, the torsionally- induced effects would 
be manifest only at the Planck scale. One can, of course, consider more general torsional 
completions of the Hilbert-Einstein action. One possibility is to introduce squared 
torsional contributions in the action so that the torsional terms will have their own 
coupling constant [T6| [T7] . In these theories torsion is present explicitly as an 
independent field and it has its own coupling constant which can be much larger than 
the Newton one, giving rise to torsional effects at much larger scales. 

In analogy with GR, one can also think to include semi-classical corrections in the 
ECSK theory so to obtain an action non linear in the scalar curvature. If we allow 
such an extension even more dramatic effects will emerge. In particular, the torsional 
contributions would acquire a coupling constant that does not only need to be fixed at 
scales that may be the nuclear one but it may also display a running behavior |18j . 

Up to now most of the analysis of these non linear ECSK theories has been focused 
on homogeneous and isotropic cosmologies, but the symmetries of these spacetimes 
greatly limit the effect of the additional torsional contributions. For this reason it 
is interesting to analyze these theories in more complicated spacetimes. One could 
think therefore to investigate the next level of generalization: the Bianchi spacetimes. 
Bianchi spacetimes are characterized by homogeneity, but present different degrees of 
anisotropies. They are classified according to the symmetry groups of the spacetime in 
nine different classes that include, for example, the Taub metric |19j. However, since the 
anisotropies are characterized by a symmetric tensor and the torsion / spin contributions 
are antisymmetric by definition, it is unlikely that their interaction will be determinant 
in Bianchi cosmologies. The one exception is constituted by the Bianchi I universes in 
which both the anisotropies and the spin terms can be represented by scalar functions 
and can therefore infiuence each other. Hence Bianchi I universes offer a relatively 
simple testing ground for the analysis of the effect of spin in spacetimes with torsion. 

In this paper we will analyze the dynamics of Bianchi I universes in the framework 
of /(i?)-gravity with torsion, in which the gravitational Lagrangian is a non linear 
function of the scalar curvature in presence of torsion. The aim is to estimate the 
influence that these non linear terms have on the behavior of the cosmology and in 
particular on the development of the anisotropies. We will perform this analysis by 
means of the Dynamical Systems Approach (DSA). This approach was proposed for the 
flrst time in its flnal form by Wainwright and Ellis and had a key role in the analysis of 
the dynamics of Bianchi universes pO] . In the last few years this technique has been 
adapted to the investigation of the cosmologies of modifled gravity revealing a series of 
interesting features [211 [221 |23l [Ml ISS [261 1271 [281 [29] . In the present paper the DSA 
will be combined with the 1+3 covarian approach. This formalism was flrst proposed 
in [301 iSll l32l [33] and it has been used to analyze a wide variety of different problems 
including the dynamics of linear cosmological perturbations [351 EEl [371 EH [391 HO] 
and the investigation of the features of modiflcations of the gravitational interaction 
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The paper will be organized as follows. In Section II we will review the general 
/(i?)-gravity theory with torsion coupled to spin fluids. In Section III we will write 
the key cosmological equations in the 1+3 covariant approach on which the dynamical 
systems approach will based. Section V will be dedicated to the phase space analysis of 
the cosmological equations. Finally, Section VI will be dedicated to the conclusions. 

Unless otherwise specified, natural units {h = c = ks = SirG = 1) will be used 
throughout this paper, Latin indices run from to 3. The symbol V represents the 
covariant derivative and d corresponds to partial differentiation in the Riemann-Cartan 
space-time [HI El HD]- We use the —,+,+,+ signature and the Riemann tensor is defined 
by 

where the T-^ are the coefficients of the linear dynamical connection, defined by 

v«.a = r."^. (2) 

The Ricci tensor is obtained by contracting the first and the third indices via the metric 

9ab 

nab = n'^acb ■ (3) 

The purely metric Riemann and Ricci tensors will be defined with the same way using 
the metric connection T^j^ and will be indicated by R^'^'^^^ , Rab, etc.. In general all the 
metric quantities will be indicated with a tilda. Finally, the symmetrization and the 
antisymmetrization over the indexes of a tensor are defined as 

T(ab) = - {Tab + Tba) , T[ab] = - {Tab " ^ba) • (4) 

2. /(7^)-gravity with torsion coupled to spin fluids 

The field equations of /(7^)-gravity with torsion are [l5l |16] 

f'{n)nab - \f{n)gab = M^b, (5) 



T„, 



f'{n) 



1 /^/'(7^) 



2 V dxP 



V) - ^a^l) + Sab' 



(6) 



In eqs. ([MS]) TZab and T^^'^ are the Ricci and Torsion tensors associated with the 
dynamical gravitational variables (fi',r), where g is a. metric tensor and F is a metric- 
compatible linear connection; Mab and 5"^^^ are the energy-momentum and spin density 
tensors of the matter fields. 
In general, the trace of equations ([5]) 

/'(7^)7^-2/(7^) = M, (7) 

is supposed to give rise to an invertible relation between the Ricci scalar curvature IZ 
and the trace M of the energy-momentum tensor. Also, it is assumed that f{TV) 7^ kTZ^ 
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(we notice that the case f{TZ) = kTZ"^ is only compatible with the condition M = 0). 
Under the stated conditions, from equation ([7]) it is possible to express the Ricci scalar 
curvature TZ suitable function of M, namely 

n = n{M). (8) 

Through the Bianchi identities |17], or also from the invariance under diffeomorphisms 
and Lorentz transformations [48], it is possible to obtain the conservation laws of the 
theory 

VaM^' + TaM'^' - MpaT'P'' ' ^5.^a7^^*"' = 0, (9) 

Moreover, we recall that any g metric-compatible connection F can be expressed as the 
sum 

r,,'^ = f,,^ - K,,\ (11) 

where K^^ are the components of the contortion tensor [8] 

KJ=\{-T,,^ + T,\-T\,). (12) 

Making use of eqs. f fTTj) and ([T2|) . from the field equations (fTOj) we obtain then the 
following representations 

KJ = k^,^ + SJ, (13) 

^ab •= ~ 2f'{7l{M)) ^^""^^ ~ ^^'^ " ~'~ ^'^"''^^ ' ^^^^ 

Kab:=-nS: + f,gP^ga,, (15) 
Tb := ..J..... fl^ + SA , (16) 



2f'{n{M)) \dxb 

which, together with equation (ITT]) , enable us to decompose the field equations (jMH]) in 
purely metric and torsional terms. 

In the present paper, we consider /(7?.)-gravity coupled to a Weyssenhoff spin fluid; the 
latter is characterized by an energy-momentum tensor of the form 

^ab ^ ^apb ^ p (^^a^b ^ ^ab^ ^ ^^^-^ 

and a spin density tensor given by 

Sab = SabU', (18) 

where {u^-Ua = — 1) and denote respectively the 4- velocity and the 4- vector density 
of energy-momentum, while Sab the spin density of the fluid (see, for example, O [10] 
and references therein). The 4- velocity and the spin density are assumed to satisfy the 
convective condition 

Sabu' = 0. (19) 
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The relations (ITU]) , together with equation ([H]), imply the identities 



0. 



(20) 



Making use of eqs. ([II]), ([I2]), ^ME^, (HE]), & and denoting by V the Levi- 

Civita covariant derivative, we can decompose the contracted curvature and the scalar 
curvature respectively as [31] 

Uab = Rab - 2Vbfa - Vcf^Qab + ^fafb - 2f,f''gab - jTcS\Ua + 

- ^ Ve {-SbaU' + S,%b - S\Ua) + S^'^ S,,UaUb, (21) 

and 

n = R- QVaf" - Qtf^ - -^^S^'^Sp,, (22) 

where now, due to the convective condition, we have Ta = ^§^- Substituting 

we get Einstein-like equations of the form [3 



expressions (12111221) into equations 



Rab - -Rgab - 



[glgt - g^'gab) VeV,/'(7^) 



2/'(7^) 



c d 

gagb 



-y^gab I X 



V J'(7^) v,/'(7^) - - (M + 7^/'(7^)) 



+ 



TqS Ma + 



+ 2f^^= + - - i77^^''^M««^'' - ^^S^'Sp^gab. (23) 

The anti-symmetrized part of the Einstein-like equations ( 123]) is equivalent to the 
conservation laws for the spin ( IIO]) 

Vc {SabU') + tS\Ua - tS^Ub + Mab - Mba = 0. (24) 

Saturating equations (IM]) with m", we get the explicit expression for the 4- vector density 
of energy-momentum 

Pb = pU,- tS\ + Ve {SabU'') U\ (25) 

where p := —u"'Pa- Inserting expressions (125]) into (fT7]) . we obtain the explicit form of 
the energy-momentum tensor 

Mab = {P + P) UaUb +pgab- UafcS",, + UaV c {SpbU") . (26) 

The symmetric part of the Einstein-like equations ( 123]) instead can be written as: 



Rab — TT-Rfi'i 



ab 



M. 



tot 
ab 



Ml + M2 + M- 



(27) 



Ml 



[glgi - /V) v,v,/'(7^) 



2/'(7^) 



glgt - i^g'^gab ] x 



VJ'(7^)v,r(7^) - 



(2J 
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= {m2 - \9a,MA , (29) 



ScdS'" -gab + UaU, , (30) 



which will be our starting point. We will define for later convenience S = ^S^-^Sab- 



3. 1+3 /(i?)-Gravity with Torsion and Spin 

In this section we will rewrite the equations of the previous section using the main tool 
of this paper: the 1+3 covariant approach. Let us start with a brief summary of the 
method. Such summary will not by any means be complete and we refer the reader to 
j30[ |3T| [32| [33] and references therein for more detailed information. 

For any given fluid 4-velocity vector field it", the projection tensor hab = gab + UaUb 
projects into the instantaneous rest-space of a comoving observer |32l [33]. The first 
covariant derivative can be decomposed as 

VaWf, = -UaUb + DaUb = - Ua Ub + I & hab + (^ab + ^ab , (31) 

where X = u'^VcX, DaXbc = h'^hlh{...V dXej..., crab is the symmetric shear tensor 

{(^ab = (^(ab), CfabU^ = 0, (T^a = 0), Uab is the VOrticity tensor {Uab = W[a6], ^^abu'' = 0) 

and iia is the acceleration vector {iia = u^VbUa). Finally, G is the volume expansion 
(O = V^Ua) which defines a length scale a along the flow lines via the standard relation 
= ^. 

a 

A general matter energy-momentum tensor Mab can always be uniquely decomposed 
with respect to u°' in the form 

Mab = IJ^UaUb + qaUb + UaQb+P hab + ^ab , (32) 

where fi is the relativistic energy density, p the isotropic pressure, the energy flux 
{qau"" = 0) and TTab the trace-free anisotropic pressure (yr^a = 0, Hab = T^iab), ^abu'^ = 0), 
all relative to u"". 

One can therefore derive the general propagation and constraint equations for the 
variables given above. These equations are called 1+3 covariant equations and their full 
set can be found in [53] . 

Here, we consider Bianchi I spacetimes. They are characterized by homogeneous 
hypersurfaces have isotropic 3-curvature Rab = \Rhab = 0. In addition, spatial 
homogeneity implies that the spatial gradients will vanish and that iia = = u. For the 
same reason qa is identically zero. Thus the characterizing kinematical 1+3 equations 
for this cosmological model are the Raychaudhuri equation 

e + |e2 + 2a2 + ^(/i + 3p) = 0, (33) 
the shear equation 

^{ab) + l^O-ab + Crc{a(Tb)'' " 77 TTab = 0, (34) 
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the definition of tlie 3-Ricci tensor 



Q = Rab= - Ol^ab) - O-ab + T^ab " \hab (2cT^ - fO^ + 2/i) , 



and its trace 



= = 2cr2 - fe^ + 2^. 



(35) 
(36) 



where Oabfy''^ = and T(^ab) = Ted 



(a b) 



Ih'^'^hab 



For our purposes it will be useful to write (135|) and ( 136|) as an equation for the shear 
propagation and a constraint between the shear, the expansion and the matter-energy 
density: 

0, (37) 
2a' + 2/i. (38) 

The first of the above equations can be used to write a dynamical equation for a: 



2 _ o„2 



a = — Bcr H . 

2 a ■ 

and will be employed in place of 
split with respect to u"- and hab'- 

+ e(/i + p) + (aV 
D^p + Db-n"^ = 0. 



The conservation equations V^'M^ 



ab 



(39) 
0, can be 

(40) 
(41) 



The set of equations (^71) . (|55]) . (HUj) form a closed system and are sufficient to 

describe completely the evolution of Bianchi I universes. 

Let us now specialize these equations in the case of /(i?)-gravity with torsion and 
spin. Like in the purely metric case [U El] this is done recasting the gravitational field 
equations in a form that resembles standard Einstein gravity plus matter and an effective 
fluid encoding the properties of the higher order part of the gravitational interaction. 
For this effective fluid the effective thermodynamics is characterized by 

52 



tot 



tot 



p 



Mlfhlu'^ 



n 



e + 



3 /"(7^)7^ 



n 



4 /'(7^) 
3 /"(7^)7^ 



3® 4 /'(7^) 



1 3(/i- + p") _ 

^ a"-^ a fin) 
/'(7^) ^ m) 



(42) 



+ 



4/'(7^) 



52 



4/'(7^)2' 



0, 



tot 



-CTab- 



(43) 
(44) 
(45) 



Using these deflnitions in the ([33]), ([37D, ([38]), ([IQ]) we obtain 



e + ie2 + 2a2 



37^ + 07^ 



//(3)(7^) r(7^)^ 



/'(7^) /'(7^)^ 
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71 3 . ^. 

■2- + 27TO^"+ ^ 



2f'{ny 



0, 



a 



a 



-fc)(7 — a— ———/<,, 



1^ f"(7Z). 
_7? _ £_i — L-j? 



e + 



3 /"(7^)7^ 



+ 



52 



3/i'^(l _ 



fi"^ + efi"'{i + w) = 0, 

5 = -65, 



9 

(46) 

(47) 



0, (4^ 



(49) 
(50) 



where the last two equations are a direct consequence of the (!9lfT0l) 

Because of the constraint ([7]) these equations can be further simphfied, but such 
simphfication can be only obtained choosing a specific form for the function /. In the 
following we will analyze the case / = 7^". 



4. TZ"- Gravity with Torsion 

We will now consider the simplest of the possible forms of / i.e. a general power of 
TZ. This model is simple enough to permit simple calculations and, at the same time 
preserves most of the most interesting features of the more general models. It will help 
us to gain an insight in the effect that torsion and spin have on the properties of the 
cosmic evolution. 

If fijVj = c/R."' then using the (171) we have 

^^/(3^^y^ (51) 



ayn — 2 

and the (HMSDI) can be written as 



Q (n — 2) n(— 2?7, + + 3)(3w — 1)" 1 An 2 

" ~ "3" (3t/;- l)(3(n- l)«; + ri-3) ^" ^ {3{n - l)w + n - 3)^ ^ 



(n - 2f-v^S''a-^''-{3w - 1)^-^^^ 
n(3(n — l)w + 71,-3) 
(nw — w — \ 



(52) 



a 



n 



Ga, (53) 



2 (3(^ — l)""^ + ^ — 3) 2 ('^ ^ 2) "tt " (3(n — l)tt) + 77, — 3)(3w — 1) " 1 

2a B H /i" + 

Gw^ n{3w — 1) 

_ (n-2)^-f«-V"(3^- 1)1-^1-^ ^2 ^ 0^ (54) 

fi"" + Qn'^{l + w) = 0, (55) 
S = -eS, (56) 

At this point we are ready to analyze the properties of the cosmology using the 
Dynamical System Approach (DSA). We will see that in spite of its original complexity 
TZ^- Gravity with Torsion can be easily analyzed with this method and we will be able 
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to highlight some interesting effect that torsion and spin generate in these cosmologies. 
Our analysis will start with the spinless case and successively we will treat the complete 
problem. 

Ji^.l. The case of zero Spin 

In the case of absence of spin the comological equations ( I52ll56l) become 

■ _ 1 2 (r2-2)-^a-^(-2Ti + 3w + 3)(3u;-l)^ i Ana'^ 

3 (3z/;-l)(3(n- 1)^ + 72-3) ^" ^ (3(n -!)«; + n - 3) ' ^ ' 

nw-w-l 

a = Bcr, (58) 

n 

[3(n-l)^ + n-3]^ 2 , [3(n - 1)^; + n - 3](3^ - 1)^ i 
Zcr — — H H J J fin = y f59j 

fi"^ + Qn"'{l + w) = 0, (60) 

In order to use the DSA, let us define the dimensionless variables 

^2 6(ri — 2)~nna~"//" (3w — 1)"~^ ^2 12n^a'^ 

~ 02(3(ri - l)«; + n-3) ' ~ e\3{n - l)w + n - 3^' ^ ' 

Note that the first of these variables is defined for 

n 3w — 1 , , 

> 0, > 0, (62) 



(3u;- l)(3(n- l)ti; + n-3) ' a{n - 2) 

if n is even or a rational number with even numerator and odd denominator; for 

n 



a{n-2){3{n- l)w + n - 3) 

if n is a rational number with odd numerator and even denominator; for 

n 



> 0, (63) 
)r and even denominator; for 

[3w — l)(3(n — l)w + n — 3) 
if n is odd or a rational number with odd numerator and odd denominator. When the 
above inequalities are not satisfied one has to redefine the first of fj6T|l with a minus 
sign. In the following we will give a unified treatment when possible. 

Using the logarithmic time = In a, the system above can be rewritten as 
jW,„f2n-3(.+ l) n^„.u3(n-l)^ + n-3^a_ 

oN [ 6n ■' 3n J 

(iS(A/') r3(n- l)«; + n-3 21 2n - 3(w + 1) 2 



dJ\f \ 3n ^ 6n ) 

1 T - = 0, (67) 

where the choice of the signs depends upon the satisfaction of the inequalities by the 
choice of n and w. Looking at the structure of the above system it is clear that the 
phase space is symmetrical for reflection and has two invariant submanifolds. The first 
is forbidden by the fact that it would correspond to absence of matter and in this case 
the (15T|) is meaningless. The second represents isotropic cosmologies. Note also that for 
n = ^^'^2"''' n = ^^i^^^ the equations decouple. 
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Using the (157|1 one can reduce this system to a single equation. If we choose the 
one for S we obtain 

^ = ^ + (1 - S^) S. (68) 

This equation admits an exact solution which shows the behavior of the "shear 
parameter" with the scale factor a 

S(a) = (69) 

.201 fjL) " +1 



where C and Sq are constant. The same can be done with the other parameter obtaining 

dVL l + W- 2nW , /r,r.\ 

W = =F ^, (1 - (70) 

which admits the solution in terms of the scale factor 

n{a) = T , ^ (71) 

It is easy to prove by direct substitution that the above solutions satisfy all of the 

In spite of these results the full dynamics of the cosmology still needs to be 
characterized by a phase space analysis. In the present case the phase space is one 
dimensional and there are two fixed points which are illustrated in Table [H 

The solutions associated with the fixed point can be found integrating the equations 
1 dQ{Af) 1 + w- nw. 



e n 
1 da{Af) _ {n-l)w~l 

a dj\f n 
1 dii{M) 



(72) 
(73) 



which are also shown in Table [1] We see, then, that A represents the shear dominated 
era and B represents an isotropic cosmology. 

Looking at the stability of these two points, A and B can be either attractors or 
repellers, but they never have the same character. Therefore our result shows that there 
are values of n for which this type of cosmology tends to dissipate anisotropics and to 
reduce to an isotropic state. Such values are given in Table |2l 

4-2. Introducing the Spin. 

We now introduce the spin of matter considering the full equations ( I52ll56l) . We could 
perform the analysis using the variables of the previous section, but it is more convenient 
to redefine them as 

6(n - 2)""?7,a~"^"(3'u; - 1)""^ _ 12n^(T^ 

" e2(3(n - l)t/; + n - 3) ' " 02(3(n - l)w; + n - 3)^ ' ^ ' 
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Table 1. The fixed points and eigenvalues for i?"-gravity witli torsion in a LRS Bianclii 
I model with spin. 

Point {Q, S) Scale Factor Energy Density Shear 

A (0,1) a = ao (t-to)='''-'"-')™J /i = a = ao{t - to)-^ 

B (1,0) a = ao {t - to)^^^ ;U = ^o(t - to) = 



Table 2. The fixed points and eigenvalues for i?"-gravity with torsion in a LRS Bianchi 
I model. 



Point (fi, S) Repeller {w = 0) Attractor(^/; = 0) 
A (0,1) n<0 n>0 
B (1,0) n>0 r2<0 

Point (fi, E) Repeller (0 < < 1) Attractor(0 < tt; < 1) 
A (0,1) n<OVn>^ 0<^<^ 
(1,0) 0<n<^ n<OVn>^ 



and naturally add the additional variable: 

^n—l<^2n—l^{n—l) 



3n 

d^^f) 3(n-l)^ + n-3 N _ ^ f]2(i-n)^i _ 2n- 3(^ + 1) ^) 

dA/" 1 3n L H J 

f/S(A/') _ ^ f 4(72 - 2) 2(2n-l)[3(n-l)^ + n-3] 2-2n 



(76) 



{3w - 1)2(^^-1) [3(n - 1)«; + n - 3]2" a202(2n-i) ' 
which accounts for the degrees of freedom of the spin. 

In this way the equations ( I52ll56l) can be written as the dynamical system 



dj\f { Qn 3n 

^ 2(3(n- l)w + n-3) ^^^2(i-n)^l 



(79) 



rfTV " 1 3 Bn 
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{2n-l)[2n-3{w + l)]^ 2(2n - l)S(3(n - + n - 3) 1 

il — > , [W) 



3n 3n 
1 - S - ^] + EQ^'^^-''^ = 0. (81 

Using the ( ISTj) . the equation for S can be ehminated 
dn{Af) 1 + w - 2nw 



dJ\f n 
dJ:{Af) 1 + w- 2nw 



(82) 

j:n, (83) 



dJ\f n 

+ - S = 0, (84) 

and it is clear that the equation for S is decoupled and can be neglected. This system 
possesses two invariant submanifolds E = and f2 = 0. For n 7^ and 1 + w — 2nw = 
the system becomes decoupled (and trivial) and can be solved exactly to give fx oc 

{n — l)w — l -. 

S oc ^ , S oc 0^ and 



which means 



= -wQ^, (85) 



(86) 
(87) 



a = 


aoit - 




t^ = 


/^o(t - 


to) , 




ao{t - 


(n — — 1 


a = 


to) 


S = 


Soit- 


'to)-'/"" 



(89) 

Using also the equation for S in which ( IHTjl is implemented 

dE(Af) l + w — 2nw,, ,^ , , 

, . / = \i2n - 1)9. - 2n + 2\E 90 

dM n ^ 

it is easy to check that the above solution satisfies all the equations. 

Setting the L.H.S of the ( I82ll83l) to zero we find the fixed points in Table [3l The 
system presents a line of fixed points and one single fixed point. The stability of these 
fixed points can be found using the Hartman-Grobman (HG) theorem and it is given 
together with their coordinates in Table |H 

The solutions for the scale factor associated to the fixed points can be found solving 
the equations 

1 dQ{U) {n-l)w-l l + w-2nw, 



Q dJ\f n 2n 

1 da{U) _ {n~l)w-l 

a dM n ' 

1 dfi{M) 



-9, (91) 

(92) 



ji dM 



1 + w), (93) 



1^ = -1 (94) 
S dM ^ ^ 

The complete solutions are shown in Table [31 

Plots of the phase space in the case of dust {w = 0) cosmologies are given in Figures 

m 
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Table 3. The solutions of the scale factor and shear evolution for i?"-gravity in a LRS 
Bianchi I model. 



Fixed Manifolds 




Scale factor 


Shear 


£ 


(0,So,0) 


n 

a = ao (t - to)^''"'" 


-iH a = ao{t-to)-^ 


A 


(1,0,0) 


2n 

a = ao (t - to)3i^+^) 


a = 




Fixed Manifolds 


(^^,S,S) 


Energy Density 


Spin Scalar 


£ 


(0,Eo,0) 


/X = 


S = So{t - to) 


A 


(1,0,0) 


li = fioit - to)-'" 


S = 



Point 
£ 



Table 4. The stability of the fixed points for i?"-gravity with torsion in a LRS Bianchi 
I cosmology with spin. Here we always assume w ^ 1/3. 



Eigenvalues 

/ r\ l+w — 2nw \ 
1^' 2n / 



l+w—2nw l+w—2nw " 



2n 



2n 



Attractor Repeller 
n>OAw=0 n<OAw=0 
< u; < 1 A (n < Vn > ^) 0<M;<lAO<n<^ 



0<n<lAO<Zi;<l 



n > 1 AO < w < 



1 

2n-l 



n<OAO<w<l 
n> 1 A <w<l 

2n—l — 




Figure 1. Phase space of the f{R) gravity with torsion in Bianchi I spacetime for 
n < 0, w — 0. The matter variable can be negative for special values of the parameters 
(see ([75])). 
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4-2.1. Asymptotics. The system f l65ll67j) is not compact therefore it is necessary to 
examine the asymptotic regime. This is done normally using the Poincare projection 
[511 [2T| . however the structure of (165116 7p allows us to use an alternative method, 
suggested for the first time in [29]. 

Starting from the (l59l) we define the variable 



D = 
and 



[3{n - l)w + n - 



|e|Vl + H2fi4(l-n)^ 



S = 2VS 



n 



(3(n — l)w + n — 3) 
6{n — 2)^^na^^{3w — 1)™ 



V ^ 33' 



2) 



(95) 

(96) 
(97) 



(3(n — l)w + n — 3) 

where the ± ensures a meaningful definition of Q for every choice of n and t/;. Defining 
the time coordinate r such that X t- = D^^X the (I65H67I) can be written as 



dr 

dQjr) 

dr 
dn\r) 

dr 
dtir) 



3w{n — 1) + n — 3^2 



11 



3n 

3w{n — 1) + n — 3 
3n 

1 + w — 2nw 



6n 



'^"+,^^-'V (i + g^) 

6n 



2r2 

{2nw — w — 1) 



2n 



(98) 

(99) 

(100) 

(101) 
(102) 



where the ± and =F are related to the definition of f2. It is clear from the previous 
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Table 5. The fixed points of the asymptotic analysis of Section 14.2.11 
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Point (Q, ^) Scale factor Shear Energy Density Spin Scalar 

~ 2n 

A- (-1,1) a = ao(t-to) ''^+™' = ;U = ;Uo(t - to)" 5 = 



2n 



^+ (1,1) a = ao{t-t^yw+^ a = fi = fio{t-to) 

£ (Qo,0) a = ao (t - to) ^^-("-d-I a = (To(t - ^o)"'^''^^''^ = 5 



system and the definitions ( I96ll97ll that if the Q is defined with a "+" sign the system 
above is compact {0 < Q < 1, < Q < 1). 

Using the constraint one can eliminate the equation for S 



(ir \ n 2n 

n\l-Q^), (104) 



gf^fi^riM, (105) 



(ir 2n 
dCl^r) 1 + w — 2nw 
dr 2n 

S2±ri2 = i. (106) 

It is clear that the first equation is decoupled so the dynamical system is only composed 
by the last three equations. The fixed manifolds can be found in the usual way and are 
indicated in Tabled One of the two fixed points {A+) have the same character of A and 
the other {A-) corresponds to the contracting versions of the solutions of A. In addition 
to A± we find an entire fixed line (£) which contains also a fixed point representing a 
power law solution that depends on the parameters n and w. The line also contains the 
origin point which corresponds to a Static Einstein universe. 

The stability of the fixed manifolds can be also calculated with the HG theorem 
and it is shown in Table ([6]). For the points common to the ones of the previous section 
the results are, naturally, consistent. The fixed line instead has a stability that depends 
on the sign of Qo('i-+w-2nw)) ^j^gj^g represents the point on the line. The point (0,0) 
has two zero eigenvalues and is not isolated. This tells us that, on the one hand, the 
point is unstable and on the other that the orientation of the orbits close to this point 
depends on the one of the neighboring points. 

Examples of the phase space are given in Figure (IHM]). 

The presence of a static solution and the fact that this solution is unstable is a sign 
that a bouncing solution is possible in this type of cosmology. A bounce in General 
Relativity is characterized by the condition 

00 = 0, Oo > 0. (107) 

holding at the bounce instant to- However in non standard gravity, since the scale factor 
does not evolve monotonically this definition is not so clearcut [52] . In fact one can only 
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Table 6. The stability of the fixed points for i?"-gravity with torsion in a LRS Bianchi 
1 cosmology with spin in the asymptotic analysis. Here we always assume w ^ 1/3. 
Note that the stability of the "+" points is the same of the one in the standard analysis 
and the one of the "-" points is reversed, as expected. 
Point Eigenvalues Attractor Repeller 

A- { -"inw+w+i -2nw+w+l ^ n< A < 10 < 1 0<n<lAO<w<l 

Ln'nJ — — — — 

1 ^ ^ 1 ^ ^ 1 A n ^ ^ 1 



> 1 A <w<l n> I < w < 



^ -in«,+u,+ M 0<n<lAO<t/;<l n<OAO<«;<l 

T^L n' nJ — — — — 

n>lAQ<w<7r-^ n > 1 A < to < 1 

— — 2n— 1 2n— 1 — 

C | (i+"'-2»^)Qo ^0| g^>o n<OAO<«;<l 0<n<lAO<u;<l 

n > 1 A -TT-^ <w <l n> I < w < 

2n— 1 — — — 2n— 1 

Qo<0 0<n<lAO<^i;<l n<OAO<u;<l 

n>lAO<ty< n > 1 A <w<l 

— — 2n— 1 2n— 1 — 




characterize if the scale factor presents a minimum, but this is not in general sufficient 
to realize the classical bounce scenario, which also require |a| <^ 1. In what follows 
we will focus only on the conditions f ll07p i.e. the necessary conditions for a bounce, 
leaving a more detailed analysis for a dedicated future work. 

From our phase space analysis it is clear that the theory presents a bouncing 
behaviour realized when the matter variable is chosen with the plus sign i.e. 

6(n — 2)~"na~" (Siw — l)^^"*^ 



(?>{n — l)w + n — 3) 



> 0, and n > 1, (108) 
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Q 

Figure 4. Phase space of the system (|103m06|) for n > 1, w = and Cl defined with 




Q 

Figure 5. Phase space of the system (jl03m06|) for n < 0, w = and defined with 



and describes a cosmology that, starting form an isotropic collapsing state bounces 
approaching an isotropic expansion state. It is interesting to note that, in a 
counterintuitive way, the matter density decreases during the collapsing phase, so that 
in the bounce it is lower than the starting value. This feature is probably due to the 
complex interaction between higher order terms and matter. Instead, when 

6(n — 2)~^na~n (^3w — 1)^^-*^ 



<0, 



and 



n < 0, 



(109) 



{3{n — l)w + n — 3) 

we have bounce connecting two shear dominated solutions i.e. a cosmology that, starting 
form an anisotropic collapsing state approaches to an anisotropic expansion state. In 
this case the matter energy density grows with the collapse, and the switch between 
expansion and contraction happens for a value of energy density higher than the initial 



one. 
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4-3. Conclusion 

In this paper we have analyzed the Bianchi I cosmology of /(/2)-gravity with torsion 
using the Dynamical Systems Approach. Compared with the standard /(i?)-gravity, 
theories with torsion contain the additional constraint ((Tj) that simplifies considerably 
the phase space analysis. In the case of f(JZ) = TU^ (E]) gives an analytical result. With 
this choice of /, we examined both the cases in which the cosmological fluid does not 
posses spin and the one in which the spin is present. 

In the first case the phase space is unidimensional and contains two fixed points 
which can be attractors or repellers. In particular, for n < and < n < '^^^ (if to 7^ 0), 
the anisotropy dominated fixed point is unstable and the matter dominated one is an 
attractor. This means that the cosmology tends to isotropize. This result is similar 
to the one found in [22] in which the case of a Bianchi I cosmology without torsion 
was analyzed. In that case, however, the limits on the values of the parameter n were 
different, but the conclusions drawn there are valid also in this case: /(i?)-gravity with 
torsion can also, provided that other observational constraints (such as the spectrum 
of gravitational waves) are satisfied, constitute an alternative to infiation. Indeed an 
isotropization effect due to torsion is already a feature of the standard EKSC theories 
of gravity, but it usually requires a contracting phase and a bounce [59] . In this respect 
the higher order corrections permit to eliminate the need for a bouncing behavior. The 
dynamical systems equations can be also solved exactly giving a useful relation between 
the matter and shear densities which could be used to perform, for example, the classical 
cosmological tests for this class of theories. 

In the case of a cosmic fluid with spin the dimension of phase space grows 
accordingly. The spin terms, however, do not enter directly in the Raychaudhuri 
equations and therefore there is no "spin dominated" flxed point. For this reason 
the introduction of the spin does not modify much the geometry of the orbits and. 
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for example, the conditions on the fixed point stabihty for the presence of an isotropic 
attractor do not change. However in this case the set of initial conditions for which this 
scenario is realized is smaller. Specifically, the isotropization of cosmic histories can only 
be realized if the initial conditions are chosen in so that < i^o < 1? which in terms of 
the standard cosmological parameters means 

1 6(n — 2)~nna~n (3w — 1)"""^ 

©V" > — — ^ — > 0- 

3(?7, — l)w + n — 3 

The cosmologies we have analyzed also include the necessary conditions for the 
occurrence of a bounce. The classical version of the bounce i.e. a phase of isotropic 
collapse followed by one of isotropic expansion is found to be possible only if n > 1 and 
the initial conditions associated with this specific behaviour are given by the inequality 
^0 > \Qq\ which, in terms of the standard cosmological quantities gives 

{n — 2)"a"(3'u; — 1) ~{?>{n — l)w + n — 3)" 



6n(3w — 1)" 
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The bouncing is also characterized by a decrease of the matter variable Vt. Since 
this variable is a monotonic function of the energy density it turns out that during 
the contraction phase the matter density actually decreases, which might appear 
counterintuitive. However one must bear in mind that there is an "effective fluid" in the 
model: the one represented by the higher order corrections. Therefore such behavior 
is most likely due to the combination of the presence of the spin and these additional 
terms. 

In addition, the phase space analysis also shows a more exotic form of bounce, 
which links a contracting and an expanding shear dominated solutions. This behavior 
is realized for any initial condition provided that 

^ — --^ ^ ^ ^ < and n < 112 

(3(n -l)w + n-?,) ^ ' 

and the matter energy density has the expected behavior i.e. grows in the contraction 
phase. It would be interesting to analyze in detail this new typology of bounce and, for 
example, understand if the anisotropy is able to "transmit" information between the 
contraction and expansion phase. 

On more general terms, our analysis conflrms, consistently with the results obtained 
in [sot [M] , that while torsion can influence in a non trivial way the behavior of f{R) 
cosmologies, the spin has a secondary role. The most likely explanation for this difference 
lays in the symmetries of the spacetime and the spin tensor, which prevent a direct 
coupling of the two quantities. If this line of reasoning is correct one can expect that 
the spin will have a key role in cosmologies in which the vorticity is non negligible, like 
Godel's universes, or, more interestingly, in the analysis of the gravitational collapse in 
which it is natural to expect an growing value of the vorticity. The exploration of these 
last cases will be an interesting further development of the present work. 
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